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Abstract 

The nonperturbative nature of nucleon-nucleon interactions evolved to low momen- 
tum has recently been investigated in free space and at finite density using Weinberg 
eigenvalues as a diagnostic. This analysis is extended here to the in-medium eigen- 
values near the Fermi surface to study pairing. For a fixed value of density and 
cutoff, the eigenvalues increase arbitrarily in magnitude close to the Fermi surface, 
signaling the pairing instability. When using normal-phase propagators, the Wein- 
berg analysis with complex energies becomes a form of stability analysis and the 
pairing gap can be estimated from the largest attractive eigenvalue. With Nambu- 
Gorkov Green's functions, the largest attractive eigenvalue goes to unity close to 
the Fermi surface, indicating the presence of bound states (Cooper pairs), and the 
corresponding eigenvector leads to the self-consistent gap function. 



1 Introduction 



The perturbativeness of a nucleon-nucleon (NN) potential can be quantified 
using the eigenvalue analysis introduced long ago by Weinberg [I]. Consider 
the operator Born series for the free-space T-matrix at energy E: 

T(E) = V + V-±—V+--- (1) 
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By finding the eigenvalues and eigenvectors of 
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V\* V )=7, V (E)\* V ), (2) 



E — H Q 

and then acting with T(E) on the eigenvectors, 

T(E)\* v ) = V\V v )(l+T h , + r£ + -~), (3) 

it follows that nonperturbative behavior at energy E is signaled by one or 
more eigenvalues with \r] u (E)\ > 1 [TH2] . Such an analysis has been recently 
used as a diagnostic of low-momentum potentials in free space [3ll^f5] . Major 
decreases in the magnitudes of the largest eigenvalues were observed as bare 
nucleon-nucleon potentials (such as those from Refs. 07]) were evolved using 
renormalization group (RG) methods [8.9]. These decreases can be linked to 
a dampening of the sources of nonperturbative physics, such as the repulsive 
core and the short-range tensor interaction, as the cutoff is lowered. Lowering 
A yields a soft potential (generically called "Viawfc")> which in turn simplifies 
few and many-body calculations [TUITT] ."^! 



The Weinberg analysis was extended to the in-medium T-matrix in Ref. [10] 
to examine the effect of Pauli blocking. The conclusion was that the bulk 
nuclear matter energy calculations are perturbative in the particle-particle 
channel. The focus in Ref. [10J was primarily on repulsive eigenvalues, which 
are tied to short-range physics. But questions naturally arise about other 
sources of non-perturbative physics. In this paper we extend the in-medium 
eigenvalue analysis to energies close to the Fermi surface, where the attractive 
NN interaction leads to a pairing instability in the particle-particle channel. 
This nonperturbative feature should be reflected in the Weinberg eigenvalues. 
Verifying the eigenvalue analysis for pairing is also a step toward a more 
general application of this tool to long-range correlations in the particle-hole 
channel. 



In Sec. [21 we review the in-medium results and focus on the Weinberg eigenval- 
ues for energies close to the Fermi surface. This takes the form of a stability 
analysis with complex energies [T3|. From the eigenvalues it is possible to 
estimate the BCS gaps. Pairing is naturally taken into account in the Nambu- 
Gorkov formalism. Using the two-particle Nambu-Gorkov Green's function 
we evaluate the Weinberg eigenvalues in Sec. [3], first for a separable model 



There are various RG methods used to generate low-momentum potentials, in- 
cluding using sharp and smooth regulators and through the Similarity Renormaliza- 
tion Group (SRG) [5J. SRG potentials depend on a parameter A, which measures the 
spread of the off-diagonal strength and acts as a cutoff. We restrict our calculations 
here to sharp-cutoff Vi ow fc potentials; gaps with smooth regulators are discussed in 
Ref. [12]. 
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and then for a Vi ow fc potential. Throughout this paper we work at zero tem- 
perature, with only a two-body potential and a free single-particle spectrum. 
Moreover we consider only the 1 So partial wave, although the analysis can be 
applied more generally and the cutoff dependence of the extracted gap used 
to test the importance of three-body forces in other partial-wave channels. We 
summarize our findings and discuss future investigations in Sec. HI 



2 In-medium Weinberg Eigenvalues and the Pairing Gap 



In free space, the Weinberg eigenvalue analysis refers to the spectrum of VGp p 
as in Eq. (j2J), where G pp is the two-body non- interacting particle-particle 
Green's function: 

As noted above, when the magnitude of the largest eigenvalue lies outside 
the unit circle, i.e., \rj u (E)\ > 1, the corresponding Born series expansion for 
the T matrix no longer converges. Lowering the cutoff results in a softening 
of the short-range repulsion and the iterated tensor force in free space, while 
the impact of shallow bound states in the S waves is eliminated at sufficiently 
high density due to Pauli-blocking. As a result the eigenvalues decrease in 
magnitude as the cutoff is lowered. In free space, for any value of the cutoff 
there are only a finite number of eigenvalues which lie outside the unit circle [T]. 

The same conclusions hold in-medium as long as we work at energies away from 
the Fermi-surface [3~][T0] . Figure [1] shows the largest attractive and repulsive 
eigenvalues at center of mass energy E = Me V for neutron matter in the 1 So 
channel for three different values of the cutoff A0 We note that lowering the 
cutoff and increasing the density both contribute to dampening the sources 
of nonperturbative physics, which is reflected by the smaller magnitudes of 
the eigenvalues at finite densities, consistent with the nuclear matter results 

of Ref. nm. 



Close to the Fermi surface, the attractive interaction between the particles 
leads to the normal ground state becoming unstable to the formation of Cooper 
pairs. While considering the behavior of the eigenvalues away from the Fermi 
surface (e.g., for analyzing bulk properties), it is sufficient to consider only 



2 At negative energy, the purely real free-space Weinberg eigenvalues can be viewed 
as inverse coupling constants that the interaction must be scaled by to support a 
bound state at that energy. Therefore, negative and positive Weinberg eigenval- 
ues are called repulsive and attractive eigenvalues. The same designations are used 
for positive-energy eigenvalues, which are complex, according to whether they are 
continuations from respulsive or attractive eigenvalues at negative energy. 
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Fig. 1. Largest attractive (rj > 0) and repulsive (rj < 0) eigenvalues at E = OMeV 
for neutron matter in the 1 S'o channel for three different cutoffs, as a function of 
Fermi momentum fer- 
tile particle-particle Green's function, because the phase space for hole prop- 
agation is small [H]. Investigating the signatures of pairing, however, requires 
us to consider the hole propagation as well, since the phase space around the 
Fermi surface becomes important. 

The in-medium Weinberg eigenvalue equation including the hole propagation 
is: 

Gl phh {E)V\S! v (E)) = Vl/ (E)\V„(E)) , (5) 

where G° pphh (E) is the in-medium two-particle and two-hole non- interacting 
propagator. In momentum space, this propagator is given by 



G° PP hh(ki, k 2 ;ui) 



N - k F )B{\k 2 \ - k F ) 6{k F - \h\)9(k F - \k 2 \) 



uj — e(ki) — e(k 2 ) + ie uj — e(ki) — e(k 2 ) 



le 



(6) 



where k\ and k 2 are the momenta of the particles (holes), k F is the Fermi 
momentum, and u is the two-particle excitation energy measured from the 
Fermi surface. The above propagator represents propagation of two particles 
above the Fermi surface and two holes below the Fermi surface. We use the 
free-particle spectrum to be consistent with the usual BCS treatment applied 
in Refs. pJ2], ^5\, and [16]. 

The phase space for pairing is maximal for back-to-back pairs [13], therefore 
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we apply Eq. ([6]) with zero center-of-mass momentum, 



0(fc - fcp) fl(fc F ~ fc) 

^pphhyk, u>) - u _ + . e - _ 2£{k) _ ^ , [ I) 



where k is the relative momentum. With uj = E — 2/x and = k 2 /2 — fi, 
where /i = kp/2 is the zero-temperature, non- interacting Fermi energy (we use 
units in which f? jm^ = 1, with the mass of a nucleon), the in- medium 
propagator is simply 

C (h F _ r.2N _ g ( fc ~ M _ g ( fc F - k ) ( R) 
^pphh\ K , 0/ 7,2 L,2 _i_ ,V P _ _ ,V ' W 



We now study the Weinberg eigenvalues for the kernel G pphh (E)V\ ow k- Just as 
for the free-space case [3], we actually solve for the eigenvalues of V\ ow kG pphh (E) , 
which has the same eigenvalue spectrum but allows for direct integration over 
singularities. In a given partial wave, the Weinberg eigenvalue equation is: 



- j\ 2 dqV M {k A ) ( j (g /J. - JP F 2 q) . ) *„(g) = . 
7r jo \ Kq — q + it Kq — q z — xe J 

(9) 

For notational convenience we have suppressed the energy dependence of the 
eigenvectors in Eq. (jSJ). Using the standard identities, 

1 1 

V Tin5(x-x ), (10) 



x — xq ± ie x — xq 



Eq. <M} becomes: 



2 ^/- A 2, T , x fO(q-k F )-0(k e -q)\ T , , 

TV JO \ /Cq — ^ y 

- zA; ^owfc(A;, fco)*„(*o) = ?7^o)^0) ■ ( 12 ) 



To identify signatures of pairing, we plot in Fig. [2] the magnitude of the largest 
attractive and repulsive Weinberg eigenvalues for G pphh V as a function of 
density for the x Sq partial wave with center-of-mass energy E = 50 MeV. 
The V[ ow fc matrix elements are from the Argonne Vi$ potential [B] using a 
sharp regulator, but the results here apply generally to all low-momentum 
potentials. As we scan through hp, the largest repulsive eigenvalue for larger 
cutoffs shows a cusp behavior due to the sharp Fermi surface that is localized 
near the momentum corresponding to kp = \^E. The repulsive eigenvalues 
are strongly cutoff dependent and no cusp is resolved at the lower cutoffs. In 
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Fig. 2. Magnitudes of the largest attractive and repulsive eigenvalues at E = 50 MeV 
as a function of Fermi momentum kp for neutron matter in the 1 Sq channel for 
three different cutoffs using the two-particle and two-hole propagator. (The repulsive 
eigenvalues are connected by lines.) 

contrast, the attractive eigenvalues depend weakly on cutoff and show a broad 
cusp near the Fermi surface; this behavior is the link to the pairing instablity. 
Similarly, Fig. [3] shows the largest attractive eigenvalue at a fixed density (in 
this case kp = 1.0 fm -1 ) as a function of E at several cutoffs. Again we find 
weak cutoff dependence and a broad cusp about the Fermi surface; in the next 
section we extract the gap from this behavior by going to complex energies. 



2.1 Stability Analysis and Pairing gaps 

In this section we adapt the stability analysis of Ref. [13]. Consider Eq. ([B]) 
for the two-particle and two-hole non-interacting Green's function. Above 2/i 
we have the particle-particle continuum and below 2/i we have the hole-hole 
continuum. Therefore, a stable bound state around 2/i cannot be accomodated. 
In fact, the bound-state energies measured from the Fermi surface 2 /a are 
purely imaginary and have a value of ±iAp, where Ap is the BCS gap at 
k-p [T3"][TT] (actually just an approximation, see below). This result can be easily 
established by studying the singularity structure of the in-medium T matrix 
in the complex plane [13]. Two purely imaginary poles of the T matrix appear 
for an attractive two-body interaction between the pairs and the normal phase 
becomes unstable. Therefore, at E = 2/i ± iA F the magnitude of the largest 
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Fig. 3. Magnitude of the largest attractive eigenvalue for neutron matter as a func- 
tion of energy E in the 1 Sq channel for three different cutoffs using the two-particle 
and two-hole propagator. 

attractive eigenvalue |r/(2/x ± iA F )\ equals one, signaling the presence of a 
bound state at these energies. Dialing the imaginary part of 2/i ± iEq from 
negative energies through to positive energies leads to eigenvalues that start 
to grow as Eq increases, cross one at Eq = ±Ap and become singular as 
Eq — > 0. The eigenvalues are symmetric about Eq = 0. 



This behavior is seen in Fig. HI where we plot the largest attractive Weinberg 
eigenvalue in the 1 5'o channel for neutron matter at k-p = 1.0 fm -1 for a cutoff 
of A = 2.0 fm -1 . The value of the imaginary part of the energy where the 
eigenvalue crosses one (E c ) directly gives a first approximation to the pairing 
gap. In Fig. [5] we show the corresponding density dependence of the pairing 
gaps extracted via the stability analysis for several cutoffs compared to self- 
consistent gaps obtained from the BCS gap equation [T2|. Note that we are 
working in the limit where Ap/fj, is small, so that the gap is independent 
of the momentum k and depends only on the density kp. The errors in the 
gaps obtained from the stability analysis scale as a power of Ap/fi. Thus the 
in-medium Weinberg eigenvalues not only reflect the instability of the normal 
phase, they also give a good estimate of the pairing gap. For the range of 
cutoffs (A = 1.6 fm~ to 2.5 fm~ ) considered here, Fig. [5] shows that the gaps 
exhibit very weak cutoff dependence, as was found in Ref. |12j . 

The gap can be cleanly extracted from the stability analysis at lower cutoffs 
because the effect of other sources of non-perturbative physics has been damp- 
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Fig. 4. Largest attractive eigenvalue as a function of energy iEq (scanning along 
the imaginary axis). E c refers to the critical energy for which \rj u \ = 1; in this case 
E c ~ Ap, the BCS gap at kp. 
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Fig. 5. Density dependence of the gap extracted through the stability analysis for 
neutron matter at several different cutoffs A compared to the self-consistent gaps 
obtained from the BCS gap equation. 
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ened, leaving the largest attractive eigenvalue isolated. The analysis is more 
involved at higher cutoffs, but still possible in the 1 So channel, where repul- 
sive eigenvalues corresponding to the strong short-range repulsion dominate. 
By continuing the attractive eigenvalue from zero energy, where it is cleanly 
distinguished, to the Fermi surface, gaps can be extracted at all cutoffs. As in 
Fig. these agree well with the BCS gaps and show only very small cutoff 
dependence at all higher cutoffs. 



3 Weinberg Eigenvalues in the Nambu-Gorkov Formalism 



Thus far we have worked in the normal phase with the non-interacting propa- 
gator. Here we consider the Weinberg eigenvalue analysis in the paired phase 
using the Nambu-Gorkov formalism. In momentum space for a homogeneous 
system, the corresponding Nambu-Gorkov Green's function reduces to [18] 

P fit y 

Gl G (k, uj) = / — \G(k, (u))g(k, + T(k, w)J*(k, uj - uj')] , (13) 

where the normal and anomalous propagators are given by 

sC^)-^^-^^—, d4) 

F{k,uj) = F ] (k,uj) = -u k v k ( — — — -) . (16) 

\uj — L k + te to + hik — le) 

The spectral functions u k and v k are defined as 



with 



E k = ^i + A(k) 2 - (18) 
Here C, k = — \i is the single-particle energy measured from the chemical 
potential /x, which for a non-interacting system at zero temperature is \i — 
kp/2. In Eq. ( TlBl . A (A;) is the gap function. Evaluating the contour integral in 
Eq. ( fi~3l) we get the following expression for the two-particle Nambu-Gorkov 
propagator, 



G° NG (k,u) = ^ V 4 . (19) 

NGV ' ' uj- 2E k + ie uj + 2E k - ie v ' 

Taking the A(k) — > limit reproduces the non-interacting two-particle hh 
propagator. 
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The Weinberg eigenvalue equation is (with E = uo + 2/i) 

G° NG (k,E)V\^(E)) = V * G (E)\*AE)) ■ (20) 

At E = 2\i the eigenvalue problem can be reduced to the BCS gap equation. 
To see this, we start with Eq. (1201) for a general potential in momentum space, 

- f A k 2 dkG ^ G (k',E)V(k',k)^ u (k,E) =rfi G (E)V v (V,E) . (21) 

TV JO 

As noted earlier, we can solve the above equation by converting it to a left- 
eigenvalue problem, or alternatively solve the right-eigenvalue problem, 

V<% Q (z)\V\9 v (z))] = »£ O (*0M**(*)>] , (22) 

which has the same spectrum as G^ G (z)V, and integrate over the singularity 
directly Using the two-particle Nambu-Gorkov propagator, Eq. fl22l) can be 
written as 



- [ A k 2 dkV(k' k) ( M fe v % \ \$ (k e) 

n Jo K ' J \u-2E k + ie ou + 2E k -ie) uK ' ; 

= r,™ G (E)V„(k',E). (23) 

For uj = (i.e., E = 2/i) and with E k = + A(k) 2 , Eq. ([23]) has a solution 
with ^ G (2/i) = 1: 

- - f A k 2 dkV{k\ k) f . 1 ] * v (k, 2/i) = 2/i) , (24) 

which is equivalent to the gap equation with the corresponding eigenvector 
ty u {k, E — > 2/i) equal to the gap function A(/c). However, it is not clear how to 
derive A(fc) from the eigenvalue equation because of the dependence of E k on 
A(/c); simple iterations do not lead to self-consistency. At energies of 2{j,±2Ap, 
the eigenvalue in Eq. fl23l) becomes singular, see Fig. [6j This coincides with 
the branch point of the two-particle Green's function and can be associated 
with the breaking of the pairs. 

To illustrate the behavior near the Fermi surface, we first calculate the Wein- 
berg eigenvalues of the Nambu-Gorkov propagator using a separable potential 
of the form 

V = \\f)(f\, (25) 
where A is a coupling that controls the strength of the potential. A is negative 
for an attractive potential and positive for a repulsive potential. The Weinberg 
eigenvalue equation in momentum space is 

-X(m u (E)) [ A k 2 dkG° NG (k,E)(k\f)\k)=^ G (E)\^ l/ (E)) . (26) 

7T JO 
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Fig. 6. Plot of the Weinberg eigenvalue for the separable potential of Eqs. (|25p 
and as a function of E for both A F = OMeV and A F k 0.604 MeV. The 
solid line represents the eigenvalue for the non-interacting propagator, G pphh , while 
the dashed lines is for the Nambu-Gorkov propagator, G^ G . The Nambu-Gorkov 
eigenvalue r]„ goes to one close to 2/i, indicating formation of a Cooper pair. Large 
singular values for rj^ G are observed at 2/i ± 2Ap, which represents the energy at 
which the pairs are broken. 

From Eq. (1261) . we see that there is only one Weinberg eigenvalue for any 
rank-one separable potential, which is given by 

V* G (E) = X(fmE)) , (27) 

and the corresponding eigenvector is 

\V V {E)) = - f A k 2 dkG° NG (k,E)f(k)\k) . (28) 

TT JO 

Substituting Eq. (128!) into Eq. (j27|) . we get the following closed expression for 
the Weinberg eigenvalue: 

V^ G (E) = A- [ A k 2 dkG° NG (k,E)\f(k)\ 2 . (29) 

TT JO 

For simplicity we choose the function f(k) to be a gaussian, 



/(*) = e~ fc2/ ' 



(30) 




Fig. 7. Largest attractive Weinberg eigenvalue as a function of energy E for a cutoff 
A = 2.0 fm -1 at kp = 1.0 fm . Notice that as E — > 2/i, r] u {E) tend to one. 

with a = \/2fm~ . We also set A = — 1 fm, so that the separable potential is 
attractive. Figure [6] shows the Weinberg eigenvalue for the separable gaussian 
potential as a function of center-of-mass energy E for Ap = OMeV and for 
finite A^. 

For vanishing Ap, the eigenvalue becomes singular at E — 2/i due to the sharp 
Fermi surface. For the gapped phase, the gap function appearing in the NG 
propagators can be estimated from the normal-phase Weinberg analysis for 
which r) v (2ji + iAp) = 1. In the present separable example, the normal-phase 
estimate gives Ap m 0.604 MeV. This value is taken as the initial guess for 
the gap function, A(k) ~ Ap, which is then used in the two-particle Nambu- 
Gorkov propagator and the corresponding eigenvalue is evaluated using 

= x\j\ W )? (^ru ~ ^fc) • (31) 



where u = E — 2/i, and E^ — J^l + Ap. We see that the eigenvalue rj^ G (E) 
tends to one close to E = 2/i, as expected. 

Figure [7] shows the largest attractive Weinberg eigenvalue as a function of E 
for V\ ow k in the 1 S partial wave at hp = 1.0 fm -1 for a cutoff A = 2.0 fm -1 . 
As with the separable example, the gap function A(k) appearing in the NG 
propagators is taken from the normal-phase Weinberg eigenvalue estimate. 
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Fig. 9. Momentum dependence of the self-consistent gap function A(fc), obtained 
from the Weinberg eigenvectors corresponding to the largest attractive eigenvalue 
at to = as the initial gap function. 
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Apart from some small numerical instabilities outside the ±2Ap region, the 
largest eigenvalue behaves like the eigenvalue in the separable E -> 2[i 

(see Fig. E]). Similarly, the eigenvector corresponding to r]^ G (E) = 1 as E — > 
2ji is an approximation to the gap function. Figure [H] shows the eigenvector 
corresponding to the largest eigenvalue at E = 2fi for representative densities. 
Using this as the first guess for the gap function A(k) and iterating the gap 
equation, 

AW= Jr A |jM, (32 ) 

yields the self-consistent BCS gap function, which is shown in Fig. [9j It is evi- 
dent that the non-self-consistent eigenvector is a poor approximation to A(k) 
in general. As kp increases, the gap closes for smaller momentum values, as 
observed in Ref. [H]. These results are consistent with the density dependence 
of the 1 Sq gap Ap. 



4 Summary 



In summary, we see that pairing instability is reflected in the behavior of the 
Weinberg eigenvalues close to the Fermi surface, with the signature of non- 
perturbativeness given by the presence of eigenvalues outside the unit circle. 
We find that a good approximation to the momentum independent gaps Ap 
can be obtained from a stability analysis, which means that weak coupling is 
a good approximation for low-momentum interactions. If instead we use the 
Nambu-Gorkov two-particle Green's function, the largest attractive eigenvalue 
tends to one close to the Fermi surface, indicating the presence of bound states 
(Cooper pairs). At the Fermi surface, the eigenvalue equation is the gap equa- 
tion and the self-consistent eigenvector corresponding to the largest attractive 
eigenvalue is the gap function. 

The gaps we have shown correspond to the BCS results and do not include 
polarization effects. These effects are known to significantly reduce the gap in 
neutron matter [T6f20lll5f21f22| . These interactions can be taken into account 
within this framework using the potential 

V(k', k) = V (k', k) + V ind (k', fc), (33) 

where V (k', k) is the two-body interaction and Vi n d(k', k) is the induced in- 
teraction, given in detail in Refs. [2"Uf2"3"] . Once again we use the two-particle 
Green's function, G° pphh , to calculate the eigenvalues close to the Fermi-surface 
and determine the in a similar fashion. This method therefore 

offers an extension to include medium effects. 

Any residual cutoff dependence in the gaps suggests the importance of many- 
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body forces. Our results agree with the recent work of Hebeler et al [12] that 
shows that the cutoff dependence of the 1 5'o gap in neutron matter is weak. 
As a result the three- nucleon contribution to the l So pairing gaps in neutron 
matter is expected to be small at the BCS level. A preliminary analysis for 
the 3 Si gaps in nuclear matter exhibit strong cutoff dependences and hence it 
would be worthwhile to investigate the role of three-nucleon interactions and 
the corresponding medium effects. 

The success of the in-medium Weinberg analysis for pairing suggests that it 
should be useful in assessing other possible sources of nonperturbative physics. 
Recent results by Roth and collaborators [23] indicate that, for low- momentum 
potentials, correlation effects in finite nuclei beyond second order do not signif- 
icantly change the binding energy per nucleon. This motivates us to investigate 
the perturbativeness of the particle-hole channel for bulk nuclear matter using 
the Weinberg eigenvalue analysis, which is in progress. 
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